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ABSTRACT 

We compare the redshifts of neutral carbon and carbon monoxide in the redshifted sources in which the 3 Pj — > 3 P fine structure tran- 
sition of neutral carbon, [Ci], has been detected, in order to measure space-time variation of the fundamental constants. Comparison 
with the CO rotational lines measures F = a 2 /p, where a is the fine structure constant and p is the electron-proton mass ratio, which 
is the same combination of constants obtained from the comparison 2 P3/2 — > 2 Pi/2 fine structure line of singly ionised carbon, [Cn], 
However, neutral carbon has the distinct advantage that it may be spatially coincident with the carbon monoxide, whereas [C n] could 
be located in the diffuse medium between molecular clouds, and so any comparison with CO could be dominated by intrinsic velocity 
differences. Using [C i], we obtain a mean variation of < AF/F >= (-3.6 ± 8.5) x 10~ 5 , over z = 2.3 - 4.1, for the eight[C i] systems, 
which degrades to (— 1.5± 1 l)x 10~ 5 , over z = 2.3 -6 A when the two [Cn] systems are included. That is, zero variation over look-back 
times of 10.8-12.8 Gyr. However, the latest optical results indicate a spatial variation in a, where Aa/a describes a dipole and we see 
the same direction in AF/F. This trend is, however, due to a single source for which the [C i] spectrum is of poor quality. This also 
applies to one of the two [Cn] spectra previously used to find a zero variation in a 2 /p. Quantifying this, we find an anti-correlation 
between \AF/F\ and the quality of the carbon detection, as measured by the spectral resolution, indicating that the typical values of 
i 50 km s , used to obtain a detection, are too coarse to reliably measure changes in the constants. From the fluxes of the known 
Z J 1 CO systems, we predict that current instruments are incapable of the sensitivities required to measure changes in the constants 
through the comparison of CO and carbon lines. We therefore discuss in detail the use of ALMA for such an undertaking and find 
that, based upon the current CO detections only, the Full Array configuration is expected to detect ~ 100 galaxies in [C i] at better than 
10 km s -1 spectral resolution, while potentially resolving the individual molecular cloud complexes at redshifts of z ~ 3. This could 
provide J 1000 individual systems with which to obtain accurate measurements of space-time variation of the constants at look-back 
times in excess of 1 1 Gyr. 

Key words. Submillimetre: galaxies - Cosmology: observations - Cosmology: theory - quasars: emission lines - galaxies: ISM - 
Techniques: spectroscopic 



1. Introduction 

Atomic carbon is an important tracer of the dense neu- 
tral gas within the Galactic and extragalactic interstellar me- 
dia. With rest freq u encies of v rest = 492.160651(55) GHz 
Yamamoto & Saitol Il99ll) and v rest = 809.34197(5) GHz 
Klein et all 119981) for the 3 P, -> 3 P and 3 P 2 -> 3 Pi 
fine structure transitions, respectively, neutral carbon is dif- 
ficult to observe at low redshift, alt hough submillime tre ob- 
servations of [Ci l in the Galaxy (|Oiha et all 1200 ll) . M82 



(IWhite et allll994IStutzki et allll997l) and other near-by galax- 
ies dGerin & PhillipsL 120001: llsrael & BaasL 120021) . show that the 



[C i] emission is closely associated with that of carbon monox- 
ide. This applies over a wide range of different environments 
and is most likely due to their similar critical densities (« ~ 
10 3 cm -3 ). Furthermore, it has been suggested that [Ci] can 
provide a measure of t he H 2 gas mass independent of CO 
ilPapadopoulos & Grevel l2004t iPapadopoulos et all 120041) . the 
standard molecular hydrogen tracer. 

Although the neutral carbon is generally believed to be lo- 
cated on the surf aces of molecular c louds, it can be found deep 
inside the cloud (llkeda et all 120021) . giving the spatial coinci- 
dences mentioned above and making the comparison of [C i] and 
CO redshifts potentially very useful in measuring the values of 
the fundamental constants at large look-back times. Comparison 



of the fine structure lines with the CO rotational lines gives 
AF/F, where F = a 2 ///, with a being the fine structure constanQ 
and yU = m e /m p the electron-proton mass ratio. Relative shifts in 
the millimetre-wave band are more than an order of magnitude 
larger than in the optical (e.g. lMurphv et al]|2001l) . from which 
the intra-comparison of redshifted ultra-violet electronic transi- 
tions, in the absorption spectra of metal-ions in damped Lyman- 
alpha system s, indicate that a has changed over the history of 
the Universe dMurphv et alll2003l) . 

As well as a temporal variation having implications for cur- 
rent grand unified theories, a spatial variation of the constants 
may address the prominent "fine tuning" issue, where if the fun- 
damental constants (which are dependent upon 27 independent 
parameters) were even slightly different than the observed val- 
ues, life could not app ear. Such a spatial variation may recently 
have been detected by I Webb et alJ d201 ll) . where a has a differ- 
ent valu es in different regions of the sky, described by a dipole 
(see also lBerengut et al1l201 ll) . 

Thus, redshifted [C i] lines are of great interest in measur- 
ing the space-time variation of the fundamental constants, espe- 
cially since the comparison CO line has already been detected 
in over a hundred cases at z ~ 1 (see Sect. 13.11 ). a number 



1 a = e 2 /he, where e is the charge on the electron, h Planck's 
constant/2^ and c the speed of light. 
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which will greatly increase with the advent of the Atacama Large 
Mi llimetre Array (ALMA) 

iLevshakov et al.l d2008l) have already compared the 2 P3/2 — » 
2 Pi/2 fine structure line of singly ionised carborQ, [C n], with the 
rotational transitions of the CO lines, which also gives F = a 2 jfi, 
at z — 4.69 and z = 6.420 They find no temporal varia- 
tion, although the uncertainties are large due to the weak [C n] 
spectra, whic h are very difficult to detect at high redshift (see 
ICurranll2009h . Furthermore, unlike the neutral carbon, ionised 
carbon traces photo-dissociation regions in the diffuse interstel- 
lar medium between molecular clouds and is therefore unlikely 
to be spatially coincident with the CO. In any case, a non- 
detected temporal variation does not necessarily imply a zero 
spatial variation, which may possibly account for discrepancies 
which arise in changes in a when attr ibuted purely to a temporal 
variation (e.g. ILevshakov et"ai1l2006l) . 

In this paper, we fit the observed [C i] and CO profiles of the 
known redshifted systems in order to derive the difference in red- 
shifts between the species in each of the sources, thus yielding 
the values of a 2 /yt at the redshifts and locations of these quasars. 



2. Analysis 

2.1. Profile fitting 

In Fig. Q] we show the 3 Pj — » 3 Po and rotational spectra of the 
sources in which both [C i] and CO have been detected^ These 
have been taken from the references listed in Table [TJ and re- 
sampled in velocity according to the reference frequency of the 
CO transition. Although the full-width half maxima (FWHMs) 
are similar between the [Ci] and CO, it is seen that in some 
cases, that there is a noticeable ve locity offset between the pro- 
files, w hich has also been noted bv lWeifi et al.l d2003l) : |Petv et al] 
d200i Pl 

We analyse the spectra in a consistent manner by fitting a 
Gaussian to each using the gaussJit. pro routine i n idlH This 
uses a Lavenburg-Marquadt algorithm dPress et al.L Il989th con- 
taining up to six terms; 



F(x) = A exp \ - 



+ A 3 + A 4 x + A 5 x 



= 1900.5369(1.3) GHz dCooksv et alii 19861) . 



3 F r om CO and [Cn] i n the quasars BR 120 2-0725 dOmont et all 
19961: llonoetall l2006bh and J 1148+5251 dBertoldi et all 12003 : 
Maiolino et alll2005l) . 

4 iLupu et all d201ll) claim a further two detections of the 3 Pi — » 3 Pi 
transition at z = 1.786 and 2.308. However, these are so poorly resolved 
as to be indistinguishable from the CO 7 — > 6 line and are thus of little 
use here. 

3 There are some differences in the FWHMs, particularly for 
H 1413+117 and APM 08279+5255, the latter of which has the poor- 
est [Ci] spectral resolution (150 km s -1 ), although the former has the 
best (17 km s~ ). The fact that in some cases the CO spectrum is wider 
than that of [Ci] and vice versa (Table [TJ, leads us to suspect that the 
differences may be due to the relatively low quality of the spectra, which 
may also affect the velocity offsets (see Sect. 12. 2\ . 

6 Although s everal rotational transi t ions have been detec t ed in some 
of the quasars {Barvainis et all 119971: iDownes et all 119991: ICox et al 
20021: iBertoldi et all 120031: IWalter et all 120031 : iTacconi et all I200L 
Ao et alll2008l : lDanielson et alll201 lh T in order to minimise differences 
in the beam filling factors, where available, we use the CO transition 
(./ = 4 — » 3) closest in frequency to the [C i] line. Except i n the case of 
H 141 3+117 (the Cloverleaf quasar), where the CO 3 -» 2 dWeifi et al 
2003) spectrum is of higher quality than the 4 — > 3 dBarvainis et al 
19971) . 



where A\ is the mean of the fit, A2 the standard deviation and 
A3, A4 and A5 fit the continuum. Since there are no measured 
uncertainties for the data, gaussJit.pro uses Poisson statistics to 
approximate data uncertainties, which is likely to overestimate 
the errors for the fits. 

Initially, fits were performed using all six terms in order 
to ensure complete removal of the continuum from the data. 
Subsequent fits were then performed, adjusting the range of the 
data being fit, as well as the number of terms used in the fit in 
order to ensure that the optimum fit was found (quantified by the 
lowest normalised^ 2 ). Generally, we found that the fits with the 
smallest errors were those with A3, A4 and A5 set to zero and, 
when applied, these coefficients had uncertainties larger than 
their fitted values. 

The results of the fits are summarised in Table Q] where the 
velocity offset (Av) is defined relative to the reference frequency 
(v re f) for each spectrum, which was then converted to a redshift 
(Zfitted)- The values for both the CO and [C 1] spectra were then 
used to derive a fractional offsetQ 

2.2. The comparison of rotational and fine structure 
transitions 

The rotational lines of CO are proportional to fi, while the fine 
structure lines of neutral or ionised carbon (represented by [C]) 
are proportional to a 2 . Comparison of the redshifts of each 
species measures 

zco - Z[C] a 2 

= , where F — — . 

zco + 1 F /i 

Showing the values for A F/F derived from t he [C i] lines, along 
with those from the [Cn] (ILevshakov e t al., 2008), we find that 
a 2 /fi exhibits little evolution with redshift (Fig. [2] left), giving 
a mean value of < AF/F >= (-3.6 + 8.5) x 10~ 5 , for the [Ci] 
systems only, and (-1.5 + 1 l)x 10~ 5 , including the [C 11] systems, 
overz = 2.3-4.1 andz = 2.3-6.4, respectively. Thus, the change 
in a 2 /\i is zero, within the observational "noise". 

Referring to Fig. [2] however, there may a decrease in AF/F 
with distance from the dipole (middle panel) and an increase 
with r cos (right panel) evident, as has been found for Aa/a 
dWebbetal.1 , 120111) . The observed range in AF/ F spans a change 
of ~ 40 times that of Aa/a (< 4 x 10~ 4 , cf . < 1 X 10~ 5 ) and 
interestingly, some theoreti cal models p redict temporal varia- 
tions of \Afi/fi\ ~ 40 |Aor/or| dUzanll2003l) . However, such mod- 
els are tentatively ruled out by recent observations, which find 
\AyM < 10~ 6 at z = 0.89 dHenkel et all l2009t iMuller et all 
[2011b and < 3 x 10~ 6 at z = 2.81 dKing etalll2011l) . 

Removing the dependence of the fine structure constant on 
the dipole [Aa/a = (1.1 + 0.2) x 10~ 6 r cos<9 - (1.9 + 0.8) x 
10~ 6 , where r is the look-back time], via, A/j/fi = 2 Aa/a - 
AF/F, gives the distribution shown in Fig. [3] However, as per 
those in Fig. [2] any apparent trend is due to the inclusion of 
APM 08279+5255 (at -10.88 Gyr and AF/F = 3.4 x 10~ 4 ), 
which has the weakest [C 1] spectrum (Fig. [T]i, while exhibit- 
ing a large offset from the CO redshift (in our fit as well as by 
Wagg et al. 2006). Indeed, several of the [C 1] detection profiles 



7 All of our analysis is done using the 3 Pi — » 3 Po transition and, 
although the 3 P 2 — » 3 Pj spectrum of [C 1] has also been detected in a 
few of the s ources - RX J09 1 1 +055 1 dWalter et alll201 lh. J11 4 8+525 1 
dAo et all 120081). SMM J1 401 1+0252 dWalter et all l20TTh. 
H 141 3+117 dWeifietall l2003h. MM 1 842 3 +593 8 dLestrade et all 
120101) SMM J2135 1 1-0102 dDanielsonetall |20TTI) . PSS 2322+1944 
dWalter et al.ll201lh . this is often blended with the CO 7 — > 6 emission. 
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Redshift, z Redshift, z Redshift, z 

2.55 2.555 2.56 2.565 2.51 2.515 2.52 2.525 3.92 3.95 3.94 




-1000 -500 500 1000 -1000 -500 500 1000 -1000 1000 

Velocity offset from z = 2.55790 [km s -1 ] Velocity offset from z = 2.51741 [km s -1 ] Velocity offset from z = 3.92960 [km s _1 ] 



Redshift, z 
2.32 2.325 2.33 



Redshift, z 
4.11 4.12 4.13 



SMM J21 351 1 -01 02 



CO 4 -» 3 
[CI] 3 P 1 -> 3 P x 1 




PSS 2322+1944 



CO 4 -* 3 

[CI] 3 P, - 3 P x 4" 




Velocity offset from z = 2.32591 [km s 



Velocity offset from z = 4.11994 [km s '] 



Fig. 1. The [Ci] 3 Pi — > 3 Po spectra (broken lines) overlaid on the CO rotational transition used (full lines) for the published high 
redshift [C i] detections, where the [C i] intensity is scaled by the value shown in the legend. All spectra are shown at the original 
spectral resolution and have been re-sampled to the reference frequency of the CO transition (v re f in TableQ] thus showing the offset 
from this), except in the case of SMMJ21 351 1-0102 (the subjec t of a spectral scan), where both spectra are shown relative to the 
derived CO 1 -> redshift of z = 2.32591 dDanielson et al.Ll201 lb. Note that for the sake of clarity, the spectra for MM 18423+5938 
have been truncated from the full +2400 km s" 1 range dLestrade et afll2010l) . 
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F = a 2 //j. 
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4 5 
CO redshift, z, 



50 100 150 

Distance from dipole, 8 [degrees] 



-10 -5 5 

Look— back time(cos 8) [Gyr] 



Fig. 2. AF/F derived from the velocity offsets in Table Q] versus the redshift (left ), distance from the dipole (middle) and the product 
of the look-back time with the cosine of the distance from the dipole (right), cf. [Webb et al. (l201lh. Th e filled black markers show 
the eight [C i] sources and the unfilled colour markers the two [C n] sources (see lLevshakov et al . 200|. 
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Table 1. Summary of the Gaussian fits to the [Ci] 3 P i — > 3 Po and the CO transitions, where we list all of the sources in 
which a [C i] detection is claimed bv lWalter et alJ (1201 ll) . although spectra are currently available for only eight (Fig.Q]). These 
are in order of increasing redshift (z quo ted, which is derived from the centroid frequencies quoted in the reference). v re f [GHz] 
is the frequency given in the reference from which the velocity offset, Av [km s ], and fitted redshift, Zatted* are determined. 
FWHM [kms -1 ] gives the full-width half maximum of the line as determined from the Gaussian fit (all single fits except 
SMM J16359+6612 for which two Gaussians were used and Av is the mean of the two fits). The final column gives AF/F, in 
units of 10~ 5 , derived from the two values of Zfuted- 





Zquoted 


Vref 


Ref 








Av 


FWHM 


Zfitted 


AF/F 


SMM J123549+6215 


2.202 


— 


T06 


CO 3 -» 


2 














— 




2.202 


— 


Wll 


[Ci] 3 P, -» 


3 Po 














— 


SDSS J102434.56 


2.2854 


105.220 


D95 


CO 3 ^ 


2 


-69 + 9 


197 


±24 


2.28565 


± 0.00010 


10 ±4 


+470909.5 


2.2854 


149.8024 


W05a 


[Ci] 3 P, -» 


3 Po 


-6.9 ± 


2.7 


157 


±31 


2.28532 


± 0.00003 




SMM J2 135-0 102 


2.32591 


138.563* 


Dll 


CO 4 


3 


+ 


13 


550 


±30 


2.32730 


± 0.00014 


-10± 10 




2.326 


147.901* 


Dll 


[Ci] 3 P, -> 


3 Po 


+ 


16 


515 


±39 


2.32764 


±0.00018 




SMM J163658+4105 


2.452 


— 


T06 


CO 3 — > 


2 














— 




2.452 


— 


Wll 


[Ci] 3 P, -» 


3 Po 














— 


SMM J16359+6612 


2.51732 


131.074 


W05 


CO 4 — > 


3 


-6 + 


21 


410+ 110 


2.51733 


± 0.00025 


-10 ± 12 




2.517 


139.9229 


Wll 


[Ci] 3 P, 


3 Po 


27 + 


15 


430 


±90 


2.51769 


±0.00018 




H 1413+117 


2.55784 


97.191 


W03 


CO 3 — > 


2 


-14 


±2 


421 ±5 


2.55774 


± 0.00002 


-4 ± 8 


(Cloverleaf) 


2.5578 


138.3313 


W05a 


[Ci] 3 Pi -» 


Po 


4 + 


22 


329 


±52 


2.55789 


± 0.00026 




SMM J1401 1+0252 


2.5653 


96.9975 


F99 


CO 3 -> 


2 


27 


±7 


189 


+ 16 


2.56532 


± 0.00008 


-14 ± 8 




2.5653 


138.0457 


W05a 


[Ci] 3 P, 


3 Po 


51 + 


18 


251 


±40 


2.56581 


±0.00021 




RXJ091 1+0551 


2.796 


91.088 


H04 


CO 3 -> 


2 


















2.796 




Wll 


[Ci] 3 P, 


3 Po 
















SMMJ02399-0136 


2.8076 


90.81 


G03 


CO 3 -> 


2 


















2.808 




Wll 


[Ci] 3 P, -» 


3 Po 
















APM 08279+5255 


3.9114 


93.867 


D99 


CO 4^ 


3 


-20 + 


12 


571 


±28 


3.91131 


± 0.00019 


-36 ±9 




3.913 


100.216 


W06 


[Ci] 3 P, -» 


3 Po 


128 + 


17 


414 


±49 


3.91309 


± 0.00027 




MM 18423+5938 


3.92960 


93.5249+ 


L10 


CO 4^ 


3 


-30 


+ 5 


162 


+ 12 


3.92911 


± 0.00009 


18 ±4 




3.930 


99.8378+ 


L10 


[Ci] 3 P, -» 


3 Po 


-83 ±6 


230 


± 14 


3.92825 


± 0.00010 




PSS 2322+1944 


4.1199 


90.048 


C02 


CO 4^ 


3 


3 + 


17 


361 


±42 


4.11999 


± 0.00029 


17 ± 13 




4.1199 


96.127 


P04 


[Ci] 3 P, -» 


3 Po 


-46 + 


22 


277 


±59 


4.11911 


± 0.00037 




SDSS Jl 148+5251 


6.4189 


93.204 


B03 


CO 6 -> 


5 


















6.419 




R09 


[Ci] 3 P 2 -> 


3 Pi 

















*The absolute frequency and +the centroid freq uency of the line from a spe ctral scan; 

References: D 95 - IDownes et aTl J1995I). D 99 - IDownes etalld 19991). F99 - iFraver et all dl999h. C 02 - ICox et alJ j2002h. G 03 - 
iGenzel et all |2003h. B03 - iBertoldi et al] |2003h. W 03 - IWeifi et al.ld2003l). H0 4 - lHainline et all d2004h. P04 - iPetv et all d2004h W05 a - 
IWeiB et al j d2005al) W05b -IWeiB et al.l d2005bl>. T06 -iTacconi et al j d2006h. W06 -IWagg et all (|2006h . R 09 -iRiechers et al J J2009h 
[according to lWalter et alj201 ill . L10 - iLestrade et al.l d2010h . Dl 1 - iDanielson et al] <201 lh . Wll - IWalter et al.U201lh . 




Look-back time(cos 8) [Gyr] 

Fig. 3. The fractional change in the electron-proton mass ratio 
versus the product of the look-back time and the cosine of the 
distance from the dipole, upon the subtracti on of 2Aa/a accord- 
ing to the best-fit dipole (We bb et"aTll201 ll) . The symbols are as 
per Fig. |2] 

are over very few channels, presumably to increase the signal- 
to-noise ratio of the spectrum. 



In order to quantify the quality of the carbon spectra, in Fig.|4] 
we show the derived value of AF/F versus the spectral resolu- 
tion, <Sv, from which we see a distinct correlation. It is surprising 
that this has a < 6% probability of occuring by chance, when ide- 
ally AF/ F should exhibit no dependence on 5v. This is strong 
evidence that any large values of AF/F in the present data are 
the result of low quality spectra. In an attempt to compensate 
for this, if we weight each value of AF/F by 1/6 v, thus having 
the lower resolution spectra making a smaller contribution to the 
mean, we obtain < AF/F > we i g hted= (-2.0 + 8.0) x 10~ 5 ([Ci] 
only) and (-1.2 + 9.3) x 10~ 5 ([Ci] and [Cn]), i.e. closer to zero 
variation. 

This interpretation of the measured values of F would 
explain why the temporal variation of F as indicated by 
APM 08279+5255 is at odds with other sources at similar red- 
shifts, which all exhibit positive values of AF/F (Fig. |2l left 
panel). As discussed above, such a discrepancy may be expected 
due to spatial variation of the constants, were the sources in 
widely separated regions of the sky. However, MM 18423+5938 
and PSS 2322+1944 are at 122 and 106 degrees from the 
dipole, respectively, whereas APM 08279+5255 is at 154 de- 
grees (Fig. [2] middle panel), thus being in the same hemisphere 
as the other two sources. This supports the suggestion that the 
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Fig. 4. The absolute value of AF/F versus the spectral resolution 
of the [C i] and [C n] detection. The symbols are as per Fig.[2]and 
the line shows the least-squares fit to all of the points. P(t) shows 
the Kendall's two-sided probability of the distribution occuring 
by chance and S (t) the corresponding significance for the [C i] 
and [Cn] spectra. 



derived offset of the [C i] redshift in APM 08279+5255 does not 
accurately trace AF/F. 

Apparent offsets may also arise from calibration errors, 
which have dominated previous measurements of the constants 
using neutral carbon: The comparison of the ultra-violet neu- 
tral carbon resonance lines with the radio-band Hi 21-cm hy- 
perfine transition yield AX/X, where X = g p a 2 /j., with g p be- 
ing the proton gyromagnetic factor. The three sources in which 
both species have been detected (all i n absorption) have yielded 
AX/X = (7 + 11) x 1(T 6 atz= 1.78 dCowie & Songailalll995l) 
and (7 + 7)xl(T 6 atz= 1.36 - 1.56 dKanekar et all 12010) . The 
relatively large uncertainties are systematic, arising from the ab- 
solute wavelength calibration of the optical spectrum. This is an 
issue which we avoid through the comparison of redshifted sub- 
millimetre transitions, since the frequency calibration of radio 
and millimetre-wave instruments is phase-locked to an atomic 
clock frequency standard. Although due to the low spectral res- 
olution of the [C i] detections and the fact that the spectra are 
emission profiles from entire galaxies, our absolute uncertain- 
ties are much larger. 

Lastly, it is entirely possible that the relative velocity shifts 
between the species are intrinsic, although as noted above, the 
values of AF/F are consistent with observational noise, with 
most of the error bars in Fig. [2] extending to AF/F » 0. 
Furthermore the FWHMs of the CO and [C i] emission profiles 
are very similar and the fact that CO is wider than [C i] in only 
about half of the cases (5 out of 8, although the uncertainties 
overlap), may indicate that the FWHMs could be the same with 
apparent differences introduced by the low spectral resolution 
of the data0 On this issue, it is also possible that systematics 
may be introduced by unresolved components in the data, again 
reinforcing the need for higher resolution data. Such data may 
also be used to ascertain the validity of fitting Gaussians to the 



The mean FWHM of the eight CO spectra is 360 (a = 150) km s" 



cf. 330 (cr ■■ 
Tabled 



110) km s for the [C i] spectra, based upon the values in 



profiles, which we believe is a reasonable assumption for the 
distribution of gas in a galaxy (with the number of G aussians 
requir ed depending on the orientation of the galaxy, e.g. lCurrar] 
12000) . However, the fact that these fits can yield "accuracies" to 
well within a single channel, gives the relatively small error bars 
for MM 18423+5938 (at z = 3.9 and AF/F = 1 .8 x 10~ 4 , Fig.0. 
Specifically, ±6 km s_ from our fi t , cf. th e ±24 km s (half a 
channel) quoted by iLestrade et all (I201Q) . Thus, an artificially 
small uncertainty could arise since, of all the [C i] spectra, this is 
the most Gaussian shaped (Fig. [TJ, when in reality the resolved 
emission may be too complex to be best fit by a single Gaussian. 
Again this is an issue which can only be addressed with higher 
resolution carbon spectra. 

To surmise, based upon the current eight high redshift [C i] 
and two [C n] detections, the data appear to be dominated by 
scatter. We have shown that this is, at least in part, due to the 
low spectral resolutions of the carbon detections. Since current 
instruments are generally capable of much finer resolutions, we 
suspect that such coarse values have been used in order to tease 
as strong a detection as possible from the data, leading to the 
scatter in AF/F. Therefore in order to make meaningful mea- 
surements of the constants, spectra with much higher signal-to- 
noise ratios are required. 



3. Improving the data 

3. 1 . Prospects with current instruments 

Of the known high redshift CO emitter^ most have the 3 Pj — > 
3 Po transition of [C i] redshifted into current millimetre bands 
and are therefore potentially detectable with current instruments, 
expanding the dataset by an order of magnitude. In order to de- 
termine the expected emission stre ngths of the [C i], from the 
current detections and those listed in W alter et aD (l201lh . we find 
a mean of I[ci] ~ .5 + 0.2 x 7co (Fig. ID, cf. the ratio s of a; 0.1 



withi n the Galaxy 



M82 dWhite et al 



Oiha et alll200lHlkeda et afll2002h . * 0.5 in 
1994tlStutzki et allll997l) "and * 0.2 in other 



near-by galaxies fcerin & PhillipsL 12000) . Since the FWHM of 
the [C i] profile is close to that of the CO in the detected cases 
(Fig. [TJ, we can replace the integrated intensities with the line 
fluxes, in order to estimate the relative telescope sensitivities 
(Fig.©. 

In the figure the expected r.m.s. noise levels achievable by 
the telescopes are derived for a spectral resolution of 50 km s 
(17 MHz at 100 GHz, typical of the current [C i] detections), af- 
ter 10 hours on source. We have used the specified system tem- 
perature of r sys = 250 K for each, which gives an r.m.s of 6 mJy 
for the Institut de Radioastronomie Millimetrique (IRAM) 30-m, 
4 mJy for the Nobeyama Radio Observatory (NRO) 45-m and 
1.3 mJy for the Australia Telescope Compact Array (ATCA). 
These sensitivities do indicate that more detections are to be ex- 
pected over all redshift ranges with current instruments. 
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Fig. 6. The flux density (S = 7/FWHM) of the CO and [C i] emitters versus redshift. The black stars show the current redshifted CO 
emitters which also exhibit [C i] emission (designated by coloured stars). The small squares show all of the other published z > 1 CO 
emitters. The telescope labels show the expected r.m.s. noise level per each 50 km s channel after 10 hours of observation with the 
corresponding t elescope at 100 GH z (see main text). The look-back time is calculated using Hq =71 km s _1 Mpc -1 , f2 ma tter = 0.27 
and — 0.73 dSpergel et al.ll2003l) . which we use throughout the paper. 
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Integrated intensity, l co [Jy km s ] 

Fig. 5. The velocity integrated intensity of the [Ci] 3 Pi — > 3 Pn 
line versus that of the CO transition quoted in TableQ] The main- 
beam temperature of MM 18423+5938 is conve rted to Jy using a 
conversion of5.0Jy K 1 dLestrade et al.Ll2010l) . The line shows 
the least-squares fit to all of the points. 

As discussed above, however, in order to measure veloc- 
ity offsets of high enough quality to measure variation in the 
constants, the spectral resolution must generally be much finer 
than for the current [C i] detections, requiring significantly lower 
r.m.s. noise levels than shown in Fig. [6] Furthermore, even with 
sufficient spectral resolution, current spatial resolutions may be 
insufficient: The finest of the current carbon detections subtend 
(unlensed) linear extents in excess of 20 kpc (HPBWs « 2", 
iDownes et al.l 1 995L 1 1 999b IWeiB et al.l2003l although > 100 kpc 
is more common. Since, we ideally require a comparison of red- 
shifts from a single complex, namely a giant molecular cloud 



(GMC), resolutions of 5 100 pc are required, which even the 
finest spatial resolutions currently available cannot resolve: For 
example, with a maximum 6-km baseline, the ATCA is capable 
of a 0.3" resolution at 100 GHz, which gives linear extents of 
between 2.2 kpc (at z = 3.7) and 1.9 kpc (at z = 4.9, the [Ci] 
redshift range observable with the ATCA 3-mm band). 

3.2. Prospects with ALMA 

On the basis of their sensitivity and spatial resolution, it is clear 
the current instruments are of limited use in using carbon to mea- 
sure cosmological variation of the constants. With the ALMA 
Early Science configuration (16 x 12-m antenna), the maximum 
400 metre baseline gives a beam of > 2" at 100 GHz, i.e. a linear 
extent of i 1 kpc at z ~ 4 and, after 10 hours on source, an r.m.s. 
noise level of 0.19 mJy per each 50 km s'channel is expected 
(Fig. [6}. For the Full Array configuration (> 50 x 12-m anten- 
nas), < 0.06 mJy per 50 km s channel will be reached after 10 
hours on source. The finest beam-width will be « 0.05 - 0.04" 
over 84-116 GHz (band-3), which corresponds to linear extent 
of 330 pc at z = 4.9 and 300 pc at z — 3.2. This level of resolu- 
tion is getting close to the S 100 pc required to resolve individual 
GMCs and, in the case where the galaxy is gravitationally lensed 
by an intervening galaxy, a magnification of only t 3 is required 
to achieve sub- 100 pc resolution. 

With ALMA, we are thus in the realm of Galactic resolution 
millimetre astronomy at very large look-back times. As such, 
the spectral resolutions used for the eight [Ci] detection ex- 
periments (Sect. 12. 2b are insufficient in the study of individual 
GMCs, as well as the application of these in measuring variation 
in the constants, where, for a ratio of line frequencies of a few 
parts in ~ 10 5 , we require spectral resolutions of better than 10 
kms" 1 . Using the ALMA sensitivity calculatoiFl at the above 
(50 km s ) and finer (10 & 1 km s ) spectral resolutions, we 
show the r.m.s. noise levels expected from the Full Array con- 



http://almascience.eso.org/call-for-proposals/sensitivity-calculator 
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Fig. 7. As per Fig.[6]but showing the expected flux density of the [Ci] emission based upon Sco ~ 2 x 5[Ci] (Fig. [5j. The dotted, 
broken and full traces show the r.m.s. sensitivity of the full ALMA array per each 1, 10 and 50 km s channel for the conditions 
used above (10 hours on source at a declination of 5 = -30° and a vapour column density of 1.262 mm). The lower redshift trace 
spans the range given by the band-7, 6, 5 and 4 receivers (370 - 125 GHz) and the higher redshift spans that given by the band-3 
and 2 receivers (1 16 - 67 GHz). The frequency scale is derived for the redshift of the [C i] 3 Pi — > 3 Po transition. The stars show the 
known high redshifted [C i] emitters with the dotted line showing the least-squares fit to these. 



figuration after 10 hours on source (Fig. [7]). From this, we see 
that the currently detected [C i] emitters are likely to be the most 
luminous of this sample, with the estimated fluxes of those yet to 
be searched generally lying below the least-squares fit described 
by these. This suggests that any further detections with current 
instruments are unlikely to be of higher quality than those in 
Fig. Q] while the ALMA Full Array will be capable of detect- 
ing [Ci] in aid of these at a spectral resolution of 10 km s or 
better. 

As seen from Fig. UJ most of the redshift space is covered, 
thus making the total number of sources expected to be de- 
tected close to the number quoted above. That is, based upon 
the current CO detections only, the Full Array configuration 
is expected to detect ~ 100 galaxies in [Ci] at better than 10 
km s spectral resolution. With regard to individual complexes, 
GMC numbers in exte rnal g alaxies can range from a few (e.g. 
6 in M31. ISheth et alj 12008) to over a hundred (148 in M33, 



lEngargiola et al Jl2003h . where observations are of sufficient sen 
sitivity. Therefore, as an order of magnitude estimate, we can 
expect the detection of J 1000 [Ci] emitters to very large look- 
back times, a much larger sample with which to measure the 
constants than available in the optical regime. 

4. Conclusions 

By fitting the carbon and carbon monoxide emission spectra, 
we have determined the relative velocity offsets between these 
species in the eight redshifted sources in which the 3 Pj — » 3 Po 
transition of neutral carbon has been detected. Such comparison 
provides a measure of the fundamental constants, specifically 
F = a 1 In, where a is the fine structure constant and fi is the 
electron-proton mass ratio. Being a comparison of fine structure 



11 Apart from six, but not on the basis of sensitivity, but that the red- 
shifted [Ci] 3 Pi — > 3 Po transitions falls between the band-4—7 and 
band-2-3 receivers at 116-125 GHz (z = 2.94 - 3.24). 



([Ci]) and rotational lines (CO), this gives the same combina- 
tion of constants as yielded by the comp arison of ionised carbon 
([C n]) with CO dLevshakovet all 120081) . although it has the dis- 
tinct advantage that the neutral carbon may be spatially coinci- 
dent with the CO, within molecular clouds, whereas the ionised 
carbon is located between clouds. 

We find < AF/F >= (-3.6 + 8.5) x 10~ 5 for the eight [Ci] 
systems and < AF/F >= (-1.5 + 11) x 10~ 5 , when the two 
[C n] systems are included, the former of which we believe to be 
the most reliable due to the segregation of [C n] and CO. Both 
results are consistent with a zero variation over the redshift range 
z = 2.3-6.4 (look-back times of 10.8-12.8 Gyr), although AF/F 
does appear to vary in the s ame direction as th e spatial dipole 
defined by the variation in a ( I Webb et al.L 1201 lb . This, however, 
is dominated by the lowest quality [C i] spectrum, which has a 
channel spacing of 150 km s , giving the largest velocity offset 
between the [C i] and CO spectra and with the removal of this 
source no trend is apparent. 

Quantifying the quality of the spectra, we find that large val- 
ues of AF/ F are favoured by the lower quality carbon spectra, 
with \AF/F\ — » as 5v — > 0, where 5v is the spectral resolu- 
tion of the carbon detection. Although the sample is small, this 
is strong evidence that the non-zero values of AF/F are the re- 
sult of low quality spectra. Conversely, however, the zero vari- 
ation of co nstants found from the ionised carbon lines are also 
unreliable dLevshakov et all 120081) by this reasoning. Weighting 
the contribution of each measurement of AF/F by \/5v, gives 
< AF/F >wei g hted= (-2.0 + 8.0) x 10~ 5 for [Ci] only and 
(-1.2 + 9.3) x 10~ 5 when the two [Cn] systems are included. 

Such low spectral resolutions are the result of maximising 
the signal-to-noise ratio in order to yield the detection. Applying 
the relation Ico ~ 2 x I[ci] to the z ~ 1 sources, indicates 
that the current [C i] detections would be among the most lu- 
minous of the sources already detected in CO. That is, even if 
current instruments were to detect [C i] emission in the remain- 
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der of the CO emitters, the spectra are unlikely to be of sufficient 
quality to provide meaningful measurements of the constants. 
Furthermore, from the spatial resolution and spectral profiles of 
the current detections, it is clear that we are detecting emission 
over much of the galaxy, when ideally we require comparison of 
the [C i] and CO lines in a single complex. As such, the best spa- 
tial resolution of 2 2 kpc, currently available in the 3-mm band 
(Zfcrj = 3.7 - 4.9), is incapable of resolving these complexes. 

The sensitivity and spectral resolution of ALMA is therefore 
required in order to fully utilise the use of neutral carbon and 
carbon monoxide in measuring the constants: 

- After 10 hours on source, the Early Science configuration 
(16 x 12-m antennas) is expected to reach an r.m.s. noise 
level of « 0.2 mJy per each 50 km s~'channel. However: 

- Although sufficient to detect [C i] in all of the currently 
known z ~ 1 CO emitters (Fig, [6] cf. Fig. |7), it is clear 
from this work that spectral resolutions of « 50 km s~ 
are insufficient. 

- Furthermore, the maximum linear resolution of i 1 kpc, 
at z ~ 4, is still too coarse to resolve individual GMCs. 

- After 10 hours on source, the Full Array configuration (> 
50xl2-m antennas) is expected to reach a similar r.m.s. noise 
level per each 10 km s -1 channel (Fig.fT). 

- This is sufficient to detect all of the current CO emitters 
in [C i] which fall into bands-3 to 7 (z = 0.33 - 6.36) at 
the spectral resolutions required. 

- The maximum angular resolution of « 0.04" at 1 16 GHz, 
gives a spatial resolution of » 300 pc at z = 3.2 - 4.9, 
which is still insufficient to resolve a $ 100 pc GMC. 
However, the majority of distant galaxies in the submil- 
limetre bands ma y be subject to magnification by grav- 
itational lensing (Blain, 19961). as a re all of the current 
rCil detections ( Ofek et al.l 120031: iGarrett et al] 120051: 
llvison etlril 12010 and references the rein) ^ Lens mag - 
nifications of > 40 may be expected dBlain et all 1 19991) . 
although magnifications of only J 3 are required for 
ALMA to achieve J 100 pc resolution at these redshifts. 

Detecting [C i] in just each of the current z ~ 1 CO emitters 
would bring the number of systems to ~ 100, which is get- 
ting close to the large sample s of optical spectra which indicate 
a spa ce-time variation in a dMurphv et all l2003t IWebb et all , 
l20ll . Being able to resolve individual GMCs in these would 
take this number to J 1000 individual systems, based upon the 
currently detected CO emitters alone. This will yield an invalu- 
able sample of high quality CO and [C 1] spectra with which 
to average out local velocity segregations which could mimic 
a change in any of the constants. 
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